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Inflation with a scalar-field potential of the form A(<ji> 2 — v 2 ) 2 can be described in terms of a 
parametrical attractor with critical points, whose driftage depends on the control value of the slowly 
changing Hubble rate. The method allows us to easily obtain theoretical expressions for fluctuations 
of inhomogeneity in both the cosmic microwave background and distribution of matter. We find the 
region for admissible values of potential parameters, wherein theoretical predictions are consistent 
with experimental results within the limits of measurement uncertainties. 
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I. INTRODUCTION 



, At present, in cosmology there is a problem in determining the parameters of inflation, which, in fact, has become 
the standard model for the early stage of Universe evolution [J, 0, H, 0, H[ before the Big Bang. The Big Bang is 
i - 1 , considered now as a short stage of reheating of the Universe due to a transformation of inflaton energy into the energy 
of matter, whereas the inflaton is usually ascribed to be a real scalar field. In this respect it would be useful to have 
a complete arsenal of effective methods in order to describe various characteristics at the inflationary stage. These 
instruments would allow us to carry out a more thorough analysis of theoretical models in comparison with quite 
precise modern experimental data. At present, the basic tool of such studies is the slow-roll approximation in the 
field equations of inflation (see the review in Q). 

The slow-roll dynamics of evolution can be consistently treated in the framework of a 1/A-expansion at a large 
amount of e-folding N for the scale factor of expansion, which was presented in Q as a general analysis of relative 
scaling behavior of inflaton quantities versus 1/N. So, the inflaton potential V gets the characteristic scale M by 
', V ~ N M 4 at M ~ 10 16 GeV, while the inflaton field <f> behaves like 4> ~ Mpi with Mpi being the Planck mass as 
t^J" ' given by the Newton gravitational constant G — l/Mp V In this respect, one could expect, for instance, a characteristic 

iy~j ■ value of quartic coupling in the inflaton self-action of the order of A ~ 1/N (M/Mp\) ~ 10~ 14 . Thus, one gets the 
<^> tool of strict consideration of the slow-rolling regime. 

0^ . We follow another way, which is the method of the quasiattractor. This approach was offered in Q for the case 
of a quadratic potential, in order to generalize and develop investigations considering the dependence of cosmological 
| evolution on initial data that were performed in [H, [i| [hJ [HLHi- Further, we have applied the same approach to the 
potential of A</> 4 in [l3[ • This kind of potential refers to the models of "chaotic inflation" , when the evolution occurs 
from large fields at Planck scales towards the global minimum at </> = 0. However, it would be useful to somehow 
generalize these results to a potential of the form \{<f> 2 — v 2 ) 2 , permitting, first, the opportunity of a situation with 
the scenario of "new inflation" , when the field evolves from a position in the vicinity of a local maximum at zero value 
of the field to the global minimum at </> — v. Second, as we will sec, such a potential allows us to essentially expand 
the region of admissible values of the potential parameters consistent with the data. This fact significantly increases 
the viability of the model. 

The quasiattractor approach can be described by the following: After choosing the model potential we derive the 
equations of the system motion, which are generally not analytically soluble, so that we try to treat the problem by 
applying some consistent approximations in order to describe the system evolution. In the method of the quasiattractor 
we introduce new dimensionless variables with presumed properties of scaling. Then, the differential equations of the 
first order can be considered as an autonomous system. The system could attain stable critical points on a phase 
plane. The trajectories converge to these points, being the attractors. Our first task is to search for such critical 
points. The notion of "quasiattractor" refers to the stable critical point of an autonomous system 1 with external 
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1 The exact attractor arises at quite definite fixed functional forms of potential: at zero cosmological constant it is the exponent [13 . H5L 
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parameters slowly drifting with the evolution. The position of the critical point is not fixed, since it is determined 
by the control parameters, which evolve and displace the critical point. But the evolution velocity of the control 
parameters is slow enough in order to consider the displacement of the point in the phase space as driftage. Thus, the 
system motion is the following: The system very quickly "falls" to the quasiattractor in the phase space, i.e. to the 
stable critical point slowly drifting during the evolution. So, the information about the initial position of the system 
is lost, while values of control parameters determining the position of attractor are important. The further evolution 
of system is exclusively determined by the driftage of the quasiattractor. The system motion is appropriated by the 
evolution of control parameters, and the system seems to lose some degrees of freedom. 

As we will derive below, the driftage of the attractor is equivalent to the slow-roll regime of inflation treated in the 
framework of 1 /TV-expansion considered in 0, . 

In Section [IT] we consider mathematical aspects of the quasiattractor, while in Section HTT1 we compare theoretical 
results with experimental data. Our results are in agreement with the precise analysis of a complete data set previously 
performed in [2(1 Hl[ in the framework of Monte Carlo Markov Chains. In the Conclusion we discuss the results 
obtained. 



II. MATHEMATICAL ASPECTS 



A. Equations of Motion 



Let us consider the action of the inflaton in the form 



s = J dx i v=g{ld„ ( f>d»<t>-v(<f>)\ (i) 

with the potential 

v = \{tf-sy. (2) 

The evolution of a homogeneous isotropic flat Universe is described by a Friedmann-Lemaitre-Robertson- Walker 
metric (FLRW) in Cartesian coordinates 

g„ v = diag(l, -a 2 {t), -a 2 (t),-a 2 {t)) 7 (3) 

where a(t) is the scale factor with its usual physical interpretation. The evolution equations read off as 

4> = -3H<j>- \<j){<j) 2 - v 2 ), (4) 
H = -4ttG0 2 . (5) 

Here the dot denotes the derivative with respect to time t. The Hubble rate is defined by H = a/a. 
The Friedmann relation is derived from ([4]) and (0), so that 

H 2 = ^^ 2 + \\{<p 2 -r?) 2 }. (6) 

All of the equations ([I])-© are consequences of General Relativity. 



B. Autonomous System 



For the sake of simplicity let us change variables: 



Il6t Il7l|. while at nonzero cosmological constant it is the hyperbolic cosine [l8j ]. 
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where k 2 = 8ttG. Then, Eqs. (g])-© take the form 2 



V\4> 2 -v 2 \, (8) 

(9) 
(10) 



I ^22 

yy = 2 X y 



3x 3 -3x- 2y 2 zy/y 2 +u 2 z 2 , (11) 
zy/y 2 + v?z 2 , (12) 



z> = ^x 2 z, (13) 

where the prime denotes the derivative with respect to N = In (a/ainit.)- Then the relation = is valid. The 

physical sense of N is that it counts the amount of e-folding during the expansion of the Universe from t- ln it. till the 
current point, i.e. when the scale factor increases by e N times. 
In terms of the new variables the Friedmann relation reads off as 

x 2 +y* = l. (14) 

The equations are simplified, since they are already differential equations of the first order, though they are nonho- 
mogeneous, which are easier for analysis than the initial ones. 

Indeed, the two equations of (fTTj) and (fT2"|) can be considered as an autonomous system of differential equations of 
the first order with external parameter z. Then, there is a question of the stability of given system. The numerical 
analysis shows, that the system is stable under some definite conditions. The control parameter of autonomous system 
is the slowly varying quantity z. Obviously, the driftage proceeds smoothly at x 2 z <C 1. The question is when will the 
critical point be stable? Then, all of trajectories will approach this point, and it becomes the parametrical attractor, 
i.e. the quasiattractor, while the system, gradually having come to it, will remain at the critical point and drift 
together with it, and the evolution of the actually stable point (the quasiattractor) will be determined by the control 
quantity z. 

Equations for the quasiattractor (x' = y' = 0) are the following: 

3x 3 -3x- 2y 2 z v / y 2 + u 2 z 2 = 0, (15) 
-xy 2 + zy/y 2 + u 2 z 2 = 0. (16) 

It is worth noticing that the system of equations is compatible with the Friedmann condition. 
For the sake of simplification of system, we can make the following change of variable: 

Y 2 = y 2 +u 2 z 2 . (17) 

Then, the equations transform and look less cumbersome without radicals: 

x' = 3x 3 - 3x- 2Yz(Y 2 - u 2 z 2 ), (18) 

Y' = ^x 2 Y + xz, (19) 

z' = \x 2 z, (20) 
1 = x 2 + (Y 2 -u 2 z 2 ) 2 . (21) 



2 For brevity of notation, we take square roots in the arithmetic sense that corresponds to the case when the field takes values greater 
than the vacuum expectation, \(j>\ > v, so that the scenario of chaotic inflation is realized. Otherwise, in the case of new inflation with 
\<j>\ < v, one should take the root with the opposite sign under the substitution y 2 — » — y' 2 in the radicand. This procedure is equivalent 
to removing the absolute value of the radicand in expression {8}, so that y 2 can formally run to negative values in the model of new 
inflation. 
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The scaling properties of Y and y are equivalent, since these quantities differ by a shift, which depends on the external 
parameter controlling the driftage. 

If the variable y is eliminated from the system, the equations for the critical point in the physical case of a; 7^ 0, 
y 7^ are reduced to the single equation in x (we recall that the quantity z is the parameter) 

^xVl - x 2 + z\j\/\ - x 2 + u 2 z 2 = 0. (22) 



C. Analysis of the System 



Let us analyze the stability of the critical point (x Cl y c ). Introduce small deviations from the attractor (Sx, 8y), then 
x = x c + Sx, y = y c + 5y. We obtain the following differential equations for the deviations in the linear approximation: 



5x' 
Sy' 



3x c y c 



\ 2 



3 Qx 2 c y A c + Au 2 z 4 

~7 ^cVc 



6y 4 c 



Sx 
Sy 



(23) 



/ 



The Friedmann condition requires 



x c 6x + 2yl5y = 0, 



(24) 



or 

( *c 2yl ) ( ^ ) = 0, (25) 
i.e. the solution, satisfying the Friedmann equation, will be proportional to the eigenvector 

Such an eigenvector for the given matrix exists, and is single with the eigenvalue 

S = 3-6 ? /-^. (27) 
Therefore, the evolution goes according to the law 

{%)= c {-iy N - < 28 » 

Thus, we see, that this is the only solution of system. It satisfies the imposed constraints. Further advancement of 
the analysis of the autonomous system will consist in the direct examination of the stability of the obtained solution. 

We require 23 < for stability of attractor. This condition is valid at small values of x and z (then, according to 
the Friedmann equation j/ 4 is close to 1) and 23 is certainly less than zero. The constraint on the smallness of x and 
z is actually valid, as we will see below. 



D. The Universe Inflation 



Let us consider Universe inflation due to the inflaton with the chosen potential. The condition of accelerated 
expansion is the following: 

a>0 => - = H + H 2 >0 => A = _3 X 2>_1. (29) 
a H z 
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Accordingly, such an expansion regime ends up with 



2 _ V3u 2 + 1-1 

Z end ~~ 



(30) 



2 1 

•^cnd g ' 

2 

vLd = 3: (31) 



(32) 



u 2 \/6 

During the actual process of expansion with acceleration, the quantities should satisfy the inequalities, 

x c < ^cndi (33) 

vi > vLi, (34) 

7? < z 2 nd . (35) 

One can see that such values are in agreement with the condition making the attractor stable (23 < 0), since 
gives 

2z 2 

2<-l--_. 

3 y 2 

hence, the accelerated expansion is governed by the stable quasiattractor. 



E. Characteristics of The Universe Expansion 

First of all, we determine how many times the Universe expands from the initial state marked by "in" to the end 
of inflation marked by "end". The total amount of e-folding TV is given by following expression, which follows from 
the equation for the parameter z in (|13[) . 

2 f z ™ d dz , s 

Nu **=*L (36) 

Here, the parameter x is set at the point of the quasiattractor. 

In order to find numerical values of the theory parameters, one should compare it with observational data. Ex- 
periment measures the inhomogeneity of the cosmic microwave background, related to the inhomogeneity of matter, 
also independently measured, hence we need to find the distribution of the inflaton inhomogeneity which leads to 
the matter inhomogeneity at the stage of reheating. Such inhomogeneity is given by the quantum fluctuations of the 
inflaton. Then, the spectral density of scalar and tensor perturbations look as 

H\ 2 fH\ 2 A 1 

and 



H\ 2 6X1 



P T (k) = 8K 2 [ — \ = — (38) 
\2ir J tx z z 4 

where the wave vector k is determined by the Hubble rate at the exit of fluctuations from the horizon, i.e. at k = aH. 
Consider the ratio r determining the relative contribution of tensor spectrum, 

r= W) = C ' (39) 

1 - (Art 

ln-*-=JV-21n — , (41) 

*^end ^cnd 

so that differentiation with respect to the wave vector is reduced to derivative with respect to the parameter z, 
determining the dynamics in the method of the quasiattractor. 



and define the spectral index ns as 



One can easily see that 
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F. Finding the Total N 



From equation ([22]) at x 2 -C 1 one approximately gets 3 



i^^(l + A 2 ). (42) 
9 

The Friedmann condition in the forms of ((6]) and (|14| strictly holds and yields y 2 « 1 in the limit under consideration, 
while the attractor position of (fT5")) . (JTHJ) reduced to (|4"2"| gives the slow- roll equation 3iT(/> + dV/dcj) — 0. 
Then, substituting the above expression into (|36p . we obtain 

3/j, 1 + u 2 ^ 2 1 

iVtotai w - u In — 

4 \ 2r 

3/4 V3u 2 +1 - 1 + V6 u 2 V6 1 2l l + u 2 zf\ 
= - u In ; ; h - U In „ „ m . (43) 

4^ v/sf+t-i \/3m 2 + 1-1 4 u 2 4 ; v ; 

To simplify this bulky expression we can separate out the function obtained by substituting z en( j, giving the term, 

F(u) = 3 - (u 2 In ^±}Z} + V ~ 6 - , ^ ) . (44) 
y ' 4^ V3u 2 + 1-1 V3u 2 + 1-1/ 

This function monotonically decreases in the interval u S [0, +oo), hence, it is restricted by limits at the borders 

~ < F{u) < \J\. (45) 

Since the inhomogeneity of the matter spectrum available for measurements actually refers to N of the order of 60, 
one can neglect the contribution of the upper limit in the integral, i.e. the value of function F(u), to the leading 
approximation in 1/iV. Then, the expression for ATtotai is simplified to 



V ^U^-uHn 1 -^^) 



Now we can express iVtotal in terms of the experimentally measured r and ns- We find 



_ 4(3z 2 +4zS 2 ) 4(94 -z 2 ) U7) 

US 4z 2 (l + u 2 z 2 )-3 3(34-1) ' 1 ' 



and express all of other parameters as follows: 



, 12r-3(r- 16) (ns - 1) 



and 



4 = 64 '- , (48) 



4 = S - (49) 



Making use of the connection between x 2 and z 2 according to (f42|) . we get 



2 64(-3r + (n s -l)(r-16)) 
u = 



(50) 



3(-4r + (n s - l)(r- 16)) 2 ' 
For the sake of simplicity, introduce the quantity % 

X = 4r-(n s -l)(r-16), (51) 



3 Let us recall, that (1421 is valid for chaotic inflation, while in the scenario of new inflation one should change the sign of y 2 , i.e. one puts 
x 2 w lz 2 (~l + u 2 z 2 ). 
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satisfying the condition of x ^ r equivalent to u ^ 0, because 



,2 



64 r — x 

Then, we can easily write down the final expression for iV to tai 



iVtotai = - 1 1 - ( 1 - - )ln (l - } . (52 ! 



XI V XJ v r ' 
The above expression can be "parametrically solved" . So, we introduce quantity ^ according to 

X = r(l-/3). 

Then 
and 

1 l-ff + fflii/? 

nS - 1 = ' iV(l - fl* -1 + /?-/? In /3 (3 + 0) - (M) 

The formula (|53|) exactly repeats the expression derived in [2(1 HH in another notation in the framework of slow-roll 
approximation, while ([54")) can match the result of 0, [2l|, if one neglects subleading terms in the denominator of 
(Ell) at l/N -> 0. 

In the limits of /3 — > and /3 — > 1 we obtain the reference cases with potentials of form X(j) 4 and m 2 <^> 2 , cor- 
respondingly. Indeed, up to corrections of the order of 1/7V as has been suggested in deriving (|55|) and ([51]) , we 
find 

hence, the vanishing of the inflaton vacuum expectation value, i.e. nullifying the quadratic term in the potential, 
while 

at tf-vl z 2 = — Ixl = — r II -31 -> 0, 
64 IAI 64 1 Ml 

that has reduced to zero the quartic term in the potential. In addition, 

f #, /? = 0, f -f, /9 = 0, 

n s -l= (55) 

I /3 = 1, I /? = !, 

in complete consistency with the consideration of these cases in other approaches. 

Generally, the scaling properties of inflation parameters versus l/N are quite complicated because of additional 
dependence on variable /3, which can correlate with the amount of e-folding N. Nevertheless, one can see that at fixed 
(J, the limit of l/N — > gives 

1 2 1 2 1 2 at 

r~ — , x ~ — , z ^ — . M~iV, 
AT' jv' 2V' 

though, actually, the dependence on /? could crucially change the asymptotic behavior: for instance, at (3 ~ expfiV] 
one gets r ~ a; 2 ~ exp[— N], and z 2 ~ u 2 ~ 1. The real situation is clarified after the appropriate analysis of the 
experimental data set. 



III. COMPARING DATA WITH EXPERIMENT 



As we have mentioned in the Introduction, a complete analysis of inflation models versus the experimental situation 
can be found in [2(| [H| , having presented the evaluation of parameters in the framework of Monte Carlo Markov 
Chains under the slow-roll approximation of theoretical entries. 



8 




FIG. 1: Data of the WMAP collaboration in the plane of the spectral parameter and the fraction of the tensor term in 
fluctuations of density: {ns,r}, in comparison with theoretical predictions at different values of e-folding N = 60 (thick solid 
line) and N — 70 (dotted line), corresponding to the exit of the fluctuation from the event horizon before the end of inflation 
(see the text). The left panel gives contours representing the WMAP data after 3 years of data taking the confidence levels equal 
to 1-cr and 2-cr, while the shaded regions give the same confidence levels after 5 years of data sampling. The right panel shows 
the WMAP data after 5 years of data taking in comparison with further constraints following from BAO and SN experiments. 



In the present paper, we will use the function in ([52")) for constructing the implicit dependence of ns versus r at 
fixed N. Theoretical curves are shown in Fig. Q]in the {ns, r}-plane (the thick solid line corresponds to N = 60, the 
dotted line does TV = 70). The dashed line shows u 2 = 0, and the region below it corresponds to the actual case of 
u 2 > 0, while the region above it marks u 2 < 0, irrelevant to the present work. 

The experimental results obtained by the WMAP collaboration after 3 and 5 years of data taking and published 
in [12] and USUI], respectively, are presented in Fig. [TJ too. The dark shaded contour gives the region with the 1-cr 
confidence level, while the shaded contour corresponds to the 2-cr level. One can see, that the theoretical calculations 
are in a good agreement with the experiment at the appropriate choice of parameters. 

From the analysis of data we can obtain quite wide limits of possible values for parameters of the model potential, 
namely 

N = 60±|g, 25 < u 2 ^ oo, (56) 

at the 1-cr level, and 

N = 60±|?, 17 u 2 sc oo, (57) 

at the 2-cr level in the 3 year data sample by WMAP. The formally infinite vacuum expectation value for the inflaton 
certainly corresponds to the final value of its mass, as we shall see below. The data acquisition of 5 year sample leads 
to more strict constraints. So, the above estimates on N with the confidence level of 1-cr transfer to the 2-cr level, as 
is clearly seen from the figure. 

However, the amount of e-folding is, in fact, limited by the actual history of the Universe evolution after inflation 
[IH, so that the analysis leads to the typical value of N m 60 (see also [Hj]). In addition, one has to take into account 
data of other experiments: that on baryonic acoustic oscillations and spacial distribution of galaxies (BAO) [26| as 
well as on the supernovae la (SN) [U HI, H^, [3(j. Such an analysis has been done in (24|, and is presented in the 
right panel of Fig. Q] 

Then, the data at the 1-cr level give the constraint on the parameter (3 in (|53"|) and (j54"| in the form 

0.75 s? sc 140. (58) 

The region of [3 ^ 1 corresponds to the scenario of chaotic inflation, when the field evolves towards the minimum of 
potential from large values at the branch of the potential approaching infinity, while (3 > 1 describes the scenario of 
new inflation, when the field "rolls down" to the minimum from small values near the peak at <fi = (one refers to 
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the case of "hilltop" inflation). Indeed, the condition for the critical point (jT6j) during inflation at y 4 — > 1 can be 
approximately written down in the form 



j ^ x 22 ^ 



at (3 < 1. So, since y 2 ~ 2 — one can straightforwardly see that /3 = 1 just separates the regions of parameters for 
new and chaotic inflation. 

Now let us determine the coupling constant A. The WMAP, BAO and SN observations give 



P S = 2.457l83g • 10- 9 , (59) 



while 



37T 2 2 

\ = 8ir 2 x 2 c z 4 P s = —r[4r-{r-16)(n s -l)] P s . (60) 
Therefore, at N w 60 with ([55)) we get 

^ A ^ 9.7 • 10~ 14 , (61) 

while the maximum is located at w 35. The scale of quartic cou plin g is quite natural, if one takes into account the 
analysis of 1/iV-expansion during the inflation as performed in [6l Il9j| and mentioned in Introduction. 
It is worth noting, that the product of Ait 2 remains finite 

Aw 2 = ^r[-3r+(n s -l)(r-16)]P s . (62) 

Moreover, under (|58[) the square of inflaton mass in vicinity of potential minimum 

2 3 2 
= 2^G AU 

takes the values 

1.03 • 10 13 GeV < m < 1.74 • 10 13 GeV, (63) 

maximal at j3 w 7. Furthermore, since sign(m 2 ) = sign(ii 2 ) = sign(u 2 ), the border of the applicability region for the 
potential is given by the following equation 

3r 

n s -l= - T „ , (64) 

16 — r 

which is represented by the dashed line in Fig. [TJ 

Experimental constraints for dependence of the spectral index on the number of e-folding N in terms of the 
parameter dns/dhxN are not restrictive, since they give a value compatible with zero at the confidence level of 2-er, 
with the quite large uncertainty being greater than the expected value of this parameter in the model under study. 
Therefore, we do not incorporate it into our estimates. 

Thus, we see, that one could extract the mass of the inflaton corresponding to maximal definiteness for all of the 
potential parameters. 



IV. CONCLUSION 



Thus, in the present paper we have carried out the analysis of an inflation model with the inflaton potential including 
both quadratic and quartic terms of self-action. The model has allowed us to consider scenarios of chaotic and new 
inflation in the framework of the quasiattractor method, which has enabled us to quite elegantly calculate the recently 
observed inhomogeneity of the cosmic microwave background and distribution of matter in the Universe. We have 
shown that such a model is consistent with the observational data. One can see, of course, that this model of the 
potential parametrically cannot satisfy all of the experimentally admissible values of ns and r within the empirical 
uncertainties (such a potential would not explain the presence of experimental points above the dashed curve in Fig. 
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[l}, but these restrictions are not critical within the accuracy of measurements, and the given potential seems to be 
consistent with the current data. 

We have obtained also, that observational data on the inhomogeneity of the Universe corresponds to the time of 
forming the inflaton fluctuations, when the Universe expands approximately e 60 times to the end of inflation, which 
is in agreement with other estimations. We have also precisely enough determined the inflaton mass. 

The work of V.V.K. is partially supported by the Russian Foundation for Basic Research, grant 07-02-00417. 
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